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Abstract 

We discuss the issue of radiation extraction in asymptotically flat 
space-times within the framework of conformal methods for numerical rel¬ 
ativity. Our aim is to show that there exists a well defined and accurate 
extraction procedure which mimics the physical measurement process. It 
operates entirely intrisically within so that there is no further approx¬ 
imation necessary apart from the basic assumption that the arena be an 
asymptotically flat space-time. We define the notion of a detector at infin¬ 
ity by idealising local observers in Minkowski space. A detailed discussion 
is presented for Maxwell fields and the generalisation to linearised and 
full gravity is performed by way of the similar structure of the asymptotic 
fields. 


1 Introduction 

In 1^] and we have described a numerical approach towards the solution of 
the Einstein vacuum equations for asymptotically flat spacetimes. It is based 
on Friedrich’s conformal field equations Q| which in turn are inspired from Pen¬ 
rose’s geometric ideas. In particular his conformal compactification procedure 
of asymptotically flat space-times allows to bring infinitely distant points in to 
finite regions. The advantages and disadvantages of this “conformal method” 
have been discussed at various places a II- Here we want to focus on one 
particular issue, namely how to extract radiation information from the numer¬ 
ically generated data and how to interpret them. This is not intended to be 
a detailed list of the routine steps to get to the relevant information because 
this has already been given in [^. Rather, we want to discuss this problem as a 
matter of principle, the emphasis of the discussion being on the mathematical 
idealizations which go into the process. 

Our motivation for doing so is mainly to clarify. People have raised some 
doubts in the past whether radiation extraction in the context of conformal 
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space-times can be a faithful process because the conformal compactification 
distorts time and space in a rather violent way thus giving rise to inaccurate 
frequencies and wavelengths of the radiation registered by a detector. Our aim 
in this paper is to present convincing arguments that, in fact, the conformal 
compactification works so beautifully that the extraction procedure based on 
the conformal method is by way of principle more faithful than procedures in 
use in standard numerical approaches. 

Certainly, most of the ideas presented in this work are not new. They are 
probably all hidden in the early works on the asymptotic structure of the zero 
rest mass fields. However, it is important to collect them here in one place 
because the issues which arise in this new conformal approach are to a large 
extent unknown to the numerical relativity community. 

The plan of the paper is as follows: we start out with a very simple radi¬ 
ating system, the Hertz dipole. We isolate the interesting quantity which is to 
be detected and show which field component is relevant. In the next section 
we discuss detectors and in particular how to idealise a detector at infinity. 
This notion is then used to demonstrate how to extract the dipole radiation on 
Minkowski space using conformal methods. Finally, we show how these ideas 
can be applied to linearized gravity and to the full non-linear theory. 


2 Radiation in Maxwell theory 

“Radiation” is a rather obscure notion in physics. Before discussing things in 
gravitational theory let us look at a simpler theory. Maxwell theory. Even there 
it is not easy to give a straightforward definition of that term. The problem is 
that “radiation” is a global concept but it is often referred to in a local context. 
One possible approach towards electromagnetic radiation is to consider a system 
as radiating if there is a nonvanishing energy flux through the sphere at infinity. 
In this statement, there is a clear reference to infinity. Another instance of the 
necessity of infinity for radiation is the fact that the radiation fields which are 
obtained from time varying multipoles are defined in the “wave zone” which in 
turn is defined as the region where the distance from the source is much larger 
than the characteristic size of the source. Hence the radiaton helds are obtained 
as a limit r oo oi the full held. In order to have an unambiguous dehnition 
of radiation one has to have access to the limit r oo. It is very important to 
state exactly how this limit is to be taken. Without specifying what happens to 
the other, especially the time coordinate, this limit may lead to unreasonable 
results. Thus, one has to study radiation from a global space-time point of of 
view. 

To get a feeling for what happens let us start out with a very simple example 
of a radiating system in Maxwell theory, the Hertz dipole. As our starting point, 
we take the relevant formulae for the helds of a radiating dipole from Jackson 0: 

B(a;,r) =/c2(n X p(w)^ — (^.l) 
r \ ikr) 

^ikr 

E(a;,r) = k‘^{n x p(a;)) x p(a;)- (2.2) 

r 

+ [3n(n.pH)-pH] 


2 RADIATION IN MAXWELL THEORY 


3 


Here, n = r/r and all hatted quantities are supposed to have a periodic time de¬ 
pendence so they are functions of the frequency io. In addition, the fields 

E and B depend on the position vector r. Unfortunately, the representation 
of the fields in terms of their Fourier transforms in time is totally inadequate 
for space-time considerations. Therefore, we undo the Fourier decomposition 
and write the fields in Minkowski coordinates (t,r). Taking into account the 
dispersion relation = 0 and choosing k = uj for retarded fields we obtain 

for the magnetic field 

B(t, r) = J B(w, r) e-**"* dto = J (n x p(w)) dto. 


Then, this expression and a similar one for the electric field result in the following 
electromagnetic dipole fields 

^ ^ P(^ “ ’’))> (2-3) 

E(t, r) = —(n X p(t — r)) X n -I —^ [3n(n • p{t — r)) — p(t — r)] 

-t ^ [3n(n • p(t - r)) - p(f - r)] . (2.4) 

We use the convention that here and in the sequel a dot will always mean the 
derivative of a function of one variable with respect to its argument. So here, p 
depends onu = t — r and the dot stands for d/du. It represents a dipole moment 
which is located at r = 0 and which varies in time. This form of the fields could 
have been obtained by applying the Lienard-Wiechert potentials to the case of 
a dipole and computing the fields from them (see e.g., Sommerfeld [l^ ). 

In order to simplify the formulae slightly, without loosing much generality we 
require p to be of the form p{u)ez, i.e., to be always pointing in a fixed direction, 
the z-axis, thus producing an axisymmetric field configuration. Introducing 
polar coordinates (r, 0, </)) and the three unit vectors Br = dr = n, eg = 1 /r ■ dg 
and e^ = l/(rsin0) • we get the simpler formulae 


E(t,r)=si„»('hd 


r2 ) 

Piu) 



eg -\-2 cos 6 




(2.5) 

( 2 . 6 ) 


This is the general field configuration for a dipole which is oriented along a fixed 
axis. It includes both the radiative fields and the field in the “near zone”. 

To extract the radiative part we consider the Poynting vector S = E x B. 
Its flux integral over a closed 2-dimensional surface represents the energy which 
flows through that surface. We obtain 

S(t, r) = sin^ 0 -I- 0(r“^) 


The exhibited term is the only one which survives integration over a sphere with 
radius R in the limit R oo. Thus, the radiated power of the system is 
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Already, we find that one has to specify what to do with the t-coordinate when 
taking the limit. Simply taking t to be constant results in an uninteresting limit 
with p being evaluated at u = — oo. We will discuss this issue in detail below. 

The radiation field consists of the terms which are proportional to p and 
these are the ones which fall off as 1/r: 

Brad(^) r) = --(n x p(u)) = sind e^, (2.7) 

1 'oiuj 

Erad(^) r) =-(n x p(u)) x n = sin6l -—- e^. (2.8) 

These are also the ones which one would obtain by solving the Maxwell equations 
for large r. By themselves they are a solution of the Maxwell equations. 

Now, what is the “radiative information” we are interested in? Obviously, 
the only information contained in the radiative fields is the function p(u). Let 
us agree, that it is this function which we are interested in, because by analysing 
the far-held, extracting p, we can obtain information about the structure of the 
source. In fact, we will show below that this is (almost) the so called “news 
function”. 

It should be emphasized that the function p here plays two roles. On the 
one hand, it describes the behaviour of the source, namely the time variation 
of the dipole, and on the other hand, it characterizes the asymptotic helds in 
the form of the news function. In the present case of Maxwell theory, both 
these properties coincide and there is a simple unique relationship between the 
far-held and the source. While this is still the case in linearized gravity it is 
no longer true in the full non-linear theory of gravity. It is expected that there 
will still be a unique relationship between the radiation helds and the structure 
(multipole moments) of the source but certainly it will not be as direct as in 
our simple example because during the way from the source to inhnity the helds 
will interact with themselves due to the non-linear nature of the theory. 


3 The detector at infinity 

The next question to consider is how to extract the information. In real life this 
is achieved by an antenna or a detector which uses the interaction with the helds 
to get information about their strength in various directions of space. Let us 
idealize such an antenna by a freely falling (non-rotating) reference frame, i.e., a 
timelike geodesic parametrized by proper time which has attached to it a triad 
of spacelike unit-vectors orthogonal to the timelike tangent and transported 
along the geodesic according to Fermi-Walker. The measuring process will be 
idealized as a multiplication of the helds with r followed by a simple projection 
onto the three axes. Then we take a detector far away from the source , i.e., 
(r, 0, (j)) hxed with r large, and try to determine the information. Since we 
cannot separate the far-held from the near-held we have to measure the full 
held at the location of the antenna and we hnd that we cannot extract a clean 
news function because there will be admixtures of the 1 /r^ and higher terms to 
the measurement result. Therefore, we need to take the antenna very far away 
from the source and, ideally, we need to take the antenna out to inhnity. 

To see more clearly what happens we now view this situation in a space-time 
which is conformal to Minkowski space, the Einstein cylinder £. To perform the 
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conformal transformation we introduce null coordinates u = t — r and v = t-\-r. 
This puts the Minkowski line-element into the form 


g = du dv — -(y — u)^da^, 


(3.1) 


where da^ is the metric of the unit-sphere. The coordinates u and v each range 
over the complete real line, subject only to the condition v — u>0. This infinite 
range is compactified by transforming with an appropriate function, e.g.. 


u = tan U, V = tan V, 


(3.2) 


thus introducing new null coordinates U and V, in terms of which the metric 
takes the form 

- TTr - TT-r{^dUdV-sin^{V-U)da‘^). (3.3) 

dcos^t/cos^F ^ ' J 

The coordinates U, V both range over the open interval (—f, f) with the re¬ 
striction V — U > 0. Obviously, the Minkowski metric becomes singular at 
points with U = ±f or V = ±f. 

Now we define a different metric g = OP'g, conformally related to g by the 
conformal factor 0 = 2 cos U cos V. Thus, 

g = MU dV-sin^iV-U)da^. (3.4) 

This metric is perfectly regular at the points mentioned above and, in fact, g is 
the metric of the Einstein cylinder. This can be verified easily by defining an 
appropriate time and radius coordinate, see |^. 

In Fig. 1^ is shown the standard conformal diagram for Minkowski space [|l^ . 
Each point in the interior of the triangle corresponds to a 2-sphere. The long 
side of the triangle consists of all the points in the center, r = 0 (i.e., U = V). 
The other two sides of the triangle correspond to the points with V = f, | ?7| < f 
{J'^) and U = —^,\V\ < ^ {U~)- These are 3-dimensional null-hypersurfaces 
which represent (future and past) “null-infinity”. The points are points in 
the center with U = V = ±f, while is a point with U = —^,V = |. These 
three points represent future and past timelike infinity and spacelike infinity. 
Thus, we may consider Minkowski space to be conformally embedded into the 
Einstein cylinder. This is shown in Fig. ^ For more information on general 
asymptotically flat space-times and their conformal extensions see [ pd| . 

The characteristic feature of the dipole held is, of course, that it is a retarded 
held: if p is such that p{u) = 0 for u < 0 then no detector will see anything at 
times t with u = t — r < 0. The information travels along the hypersurfaces 
u = const., the future light cones emanating from the world-line of the source. 
These are shown in hg. |l| as straight lines reaching out to . 

What happens now if we take a detector and move it further and further 
out to inhnity in order to extract the news function as cleanly as possible? This 
is not as straightforward as one would expect. The world-line of a detector is 
given by (ignoring the angular variables) 

7 p(A) = (t(A),r(A)) = (T + A,i?). 

The event P = 7p(0) which characterizes the world-line of the detector uniquely 
can be interpreted as the event when the detector “switches on”, i.e., starts to 
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Figure 1: The conformal diagram of Minkowski space. The lines meeting at 
are lines of constant t, while the lines emanating from i~ and converging into 
are lines of constant r. They are the world-lines of detectors at a finite distance 
from the source. The lines at 45° symbolise radiation. 



Figure 2: The embedding of Minkowski space into the Einstein cylinder £. (This 
figure is taken from 0) 
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register. We have T = t{P) and R = r{P). The world-line of the detector in 
terms of (u, r;)-coordinates is 

7p : (m(A), v{X)) = {T-R + X,T + R + X) (3.5) 

and in terms of {U, y)-coordinates one has 

7 P : {U{X),V{X)) = (arctan(T — R + X), arctan(T -|- i? -|- A)) (3.6) 

We take a sequence of detectors, determined by their “switch-on events” Pn 
(with Tn = t{Pn) and i?„ = r{Pn)) such that Rn oo for n —> oo. First 
suppose the sequence is such that Tn = T = const. Then we obtain a “limiting 
detector” 

lim 7 p„(A) = lim (C/„(A), K(A)) = (-^,^), (3.7) 

n—>-oo n^oo Z A 

i.e., the limiting detector has been compressed entirely to the point This 
limit is reasonable but not useful. In fact, what happens is that the detectors 
switch on at the same time t = T. Then they have to sit and wait for the signal 
to arrive, the delay growing longer and longer, depending on their distance to 
the source, and in the limit there will be no signal at all. 

Therefore, we need to use a different sequence of detectors. Operationally 
speaking, we trigger the detectors to switch on when the signal arrives. Thus, 
we have a sequence for which T„ — = uq = const. Then we obtain 

lim 7 p„(A) = lim (C/„(A), 14(A)) = (arctan(uo -I- A), ^). (3.8) 

n—*-cx) n —>oo 2 

This “limiting detector” is a line on determined by the angles 6 and </>. 
Furthermore, it is a generator of and, hence, a null geodesic. However, the 
parameter A is not just any parameter along that geodesic. Instead, it is a so 
called Bondi parameter |Tl|] : let s be an arbitrary parameter along the generator 
and let k = d/ds be the tangent vector along the generator with k(s) = 1. The 
acceleration a of k is defined by the equation dk/ds = ak. Then A is a Bondi 
parameter if it satishes the equation 

Here, p' is the expansion of 4^, a quantity which determines how quickly neigh¬ 
bouring generators spread apart. It is determined intrinsically from the way 
4^ is attached to the physical space-time. The concept of a Bondi parameter is 
conformally invariant so it does not depend on which conformal factor has been 
used to attach the boundary to physical space. If A is a Bondi parameter then 
every other Bondi parameter has the form G{6, (/))A -I- H{d, </>) where G and H 
are arbitrary functions on the sphere of generators of 4^ ■ 

This freedom in the choice of parameters is exactly what one would expect 
considering the way it was defined. Recall that we took the sequence of detectors 
in such a way that their world-lines were defined by keeping the spatial coordi- 
nartes fixed. However, these coordinates were one special choice of Minkowski 
coordinates. Any other choice would serve the same purpose. Thus, we have the 
freedom of the full Poincare group available to define the detectors. The limit 
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of these detectors will always be a generator of with a Bondi parameter on 
it but the parameters will not be the same. Instead they will exhaust the full 
2-dimensional affine space of possible Bondi parameters. 

To verify that A is indeed a Bondi parameter in our case we observe that 
U can be taken as a parameter along the generator. In fact, it is an affine 
parameter so that the corresponding acceleration van ishes. For the expansion we 
get p' = tan 17. Inserting this into the equation (^) shows that it is identically 
satisfied. In fact, we may take the ret arded time u = tanC/ itself as a Bondi 
parameter because it also satisfies (^). In order to complete the notion of a 
detector at infinity we observe that the radiation fields are both tangent to the 
spheres of constant t and r and this remains so even in the limit r —>■ oo. Thus, 
we require that the detector at infinity has two of its axes tangent to 77“*' and 
perpendicular to the generator. 

What we have found with this limiting procedure is a concept which is 
entirely intrinsic to 77^. It does not depend anymore on the fact that we have 
started with Minkowski space. Instead it is a concept which can be introduced 
whenever the notion of null-infinity is well defined. Therefore, we now make the 


Definition 1 A detector at infinity is idealized by a generator of which 
is parametrized by a Bondi parameter u. The axes of the detector should be 
oriented so that one axis (the “retarded time” axis) coincides with d/du while 
two other axis are perpendicular to the generator and tangent to . 


4 Radiation extraction at infinity 

Let us see whether we can now use this definition to extract radiative informa¬ 
tion. Suppose we were to solve the Maxwell equations numerically in Minkowski 
space. And suppose we do not use the traditional way but instead solve a confor¬ 
mally equivalent set of equations on a conformally related space like the Einstein 
cylinder E. The Maxwell equations are conformally invariant so we still have to 
solve the Maxwell equations, only reexpressed in terms of the conformal metric 
g instead of the Minkowski metric g. Of course, we have to put in somewhere 
that we are really interested in fields on Minkowski space and not in the more 
general fields on the Einstein cylinder £. This is done by specifying the initial 
data appropriately. 

The numerical solution of these equations provides us with a Maxwell field 
on the complete Einstein cylinder including that piece which corresponds to 
Minkowski space and its conformal boundary 77+. Usually, this solution will be 
represented on 3-dimensional spacelike hypersurfaces which form a foliation of 
the conformal space. On each such slice, we can now search for the location of 
its intersection with 77+, a 2-dimensional closed surface. Then we can evaluate 
the Maxwell field on that surface and so, keeping the angles fixed and going 
from one slice to the next, we obtain the value of the Maxwell field along a 
generator of 77+ as a function of the time parameter which labels the slices. 

As an example, we consider again the dipole fields. In order to get the 
solutions onto the Einstein cylinder we first write the electric and magnetic 
field as one-forms and then combine them into the Earaday 2-form according to 


F = dt A E + irB. 


( 4 . 1 ) 
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The star is the Hodge star of the metric on the euclidean spaces t = const. The 
covariant form of the fields turns out to be 

E = sin 9d9 + 2 cos (^ (^-^ + (4-2) 

B= sin^ 9 d(j). (4.3) 

In terms of the double null coordinates u and v we obtain 

F = — ( p + -] du Ady — {du + dv) A dy + y ( ^ + ^] du Adv, (4.4) 
\ rj 2C CJ 

where r = (v — u)/2 is considered as a function of u and v and where we have 
introduced y = cos 9 as another coordinate. Observe, that the radiation field sits 
entirely in the duAdy component of F. It is the only component which survives 
the limit r —> oo, keeping u fixed. In the non-axisymmetric case we would 
also get a radiation component proportional to du A d(j) and this is the general 
behaviour for Maxwell fields on asymptotically flat space-times. This implies 
that we get the radiation fields from the Faraday 2-form simply by restricting 
it to 

The Faraday 2-form is conformally invariant. This implies that one gets 
the corresponding Faraday form on the Einstein cylinder by performing the 
coordinate transformation which transforms from the u, v coordinates to the 
[/, V coordinates. This results in 


F = - 


— 2 tan I/P^ cos^ U + 2 
P 


P cos U cos V 


2sin^(E- [/) 


sin(E — U) 
(cos^ VdU + cos^ UdV) A dy 


dU A dy 


4y 


2PcosUcosV Pcos^P \ 

-rr— + n„. -rr- 1 dU A dV, 


sirTiV-U) sm\V-U) 


where the function P is defined by P{U) = p(tan C7) and the dot here means 
d/dU. 

Let us now restrict this form to by pulling it back along the inclusion 
map {U,y,(j)) ^ {U, ^,y,^) 

i*F = — ^P — 2 tan UP^ cos^ U dU A dy. (4.5) 

According to the discussion above, this is the radiation field. 

The detector at infinity is a line y = const, (and (j) = const.) on with 
a Bondi parameter on it. The variable U is not a Bondi parameter. In order 
to find one we need to solve the second order equation (^). Of course, we 
know that u = tan P is a solution, so that the radiation field in terms of this 
parameter reads 


i*F = —p{u) du A dy, 


(4.6) 


where P{U) = p(tanP). By comparison with the above definition of P we get 
p{u) = p{u). Hence, the radiation field obtained from the conformal treatment 
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by restriction to and introduction of a (special) Bondi parameter exactly 
reproduces the radiative information that we are interested in. 

However, in most cases of interest there is no preferred Bondi parameter and 
then we have to regard all parameters as equivalent. To see what changes when 
a different parameter is used, we recall that any Bondi parameter is of the form 
u = Gu + H. Since the shift of origin is of no consequence we choose H = 0. 
Then u = u/G and 


i*F =p{u/G) du A dy, (4.7) 

Ct 

so that the radiation field expressed in terms of u is ■^p{u/G). 

In order to interpret this behaviour we have to keep in mind that the elec¬ 
tromagnetic field has been given with respect to one specific Minkowski coor¬ 
dinate system. But all such systems are equivalent. Suppose we choose a new 
Minkowski coordinate system which is boosted with velocity w with respect to 
the old one along the direction towards the detector. Then one finds the trans¬ 
formation of the retarded time to he u u = au with a = 1^(1 -I- ic)/)! — w). 
The radiative components of the covariant fields transform according to 

('^rad’'^rad) (-^rad’^^rad) ~ ~(^rad’-®rad)- (4-8) 

For we then have 

^rad “^rad = dy = --p{u/a) dy (4.9) 

and similarly for Identifying G and a we have the same transformation 

law. This means that the freedom in the Bondi parameters reflects the freedom 
in the choice of Minkowski coordinates. Mathematically, this is to be expected 
and physically, this is a well known phenomenon, namely the Doppler shift in 
the incident signal due to the relative motion of the source and the detector. 
The ambiguity arises because this motion is not known a priori but has to be 
determined by other means (e.g., spectroscopic methods). 


5 Gravitational radiation 

Let us now see how the above discussion applies to the more complicated case 
of gravity. We have already mentioned above that the concept of the detector 
at infinity as stated in section ^ is entirely intrinsic to 77. It can be applied 
whenever there is a well-defined notion of null-infinity available. Therefore, we 
only need to see what corresponds to the radiative information which, in our 
simple example of the Maxwell case, was given by the changing dipole moment. 
For the following discussion we refer to the work of Janis and Newman on 
the structure of the sources of the gravitational field. This paper contains a 
unified treatment of Maxwell theory and linearized gravity which is very useful 
for our purposes. Incidentally, these authors treat the radiation of the general 
multipole moments from which it becomes obvious that our discussion which 
was restricted to the dipole case generalizes in a straightforward way without 
any change in the essential statements. 
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We put the electromagnetic dipole field in a form which is very much like 
the one for the gravitational field. We introduce the Bondi frame of (complex) 
null vector fields 

l=2dy, n = du, m = ^-(dg + (5.1) 

V2r \ sine 

Then we can encode the electric and magnetic fields in three complex scalar 
functions 4>2 as follows: 

())o = F(l,m) =sin6», (5.2) 

<(>i = ^ (^(1. n) + F(m, m)) = - cos 9, (5.3) 

<)>2 = F(m,n)=-^(^^ + ^ + ^^sin0 (5.4) 

A simple definition of oi = p/v^ achieves complete agreement with the dipole 
field given in j|]. These authors define the coefficient 4>%(u,6,4>) of 1/r in the 
scalar (j )2 as the “news function”. Physically, this function is the information 
sent by a “broadcasting” station. Mathematically, it is one piece of the free 
data for the solution of the characteristic initial value problem for the Maxwell 
equations on Minkowski space. It is that piece of data which can be considered 
as being given on or, physically, that component of the field which can be 
extracted from ^7+ by an infinite detector. We will also call it the “null-datum” 
on J+. 

The similarity between the Maxwell equations and the equations for lin¬ 
earized gravity suggests that also in the latter case there is a “news function”. 
Indeed, the function in question is the component of the linearized Weyl 
spinor which is entirely intrinsic to It plays the same role for the solu¬ 

tion of the characteristic initial value problem as does the ^2 the Maxwell 
case. Thus, in the linearized gravity case things are very much in analogy to 
the Maxwell case. In order to extract the radiative information one determines 
the null-datum '0° along all the generators of J''*' as a function of a Bondi 
parameter. 

For the full theory we can again employ an analogy. The asymptotic struc¬ 
ture of the non-linear gravitational field is very similar to that of the asymp¬ 
totic linearized field. Again, one can isolate a function 4^°, the component of 
the Weyl spinor which is intrinsic to This can also be given freely in the 
solution of the asymptotic characteristic initial value problem (see e.g., |^). 
Thus, mathematically, it serves the same purpose as 0° '04- However, its 

physical meaning is not so clear cut as before, dt® does give information about 
the multipole structure of the gravitational radiation but the implications for 
the multipole structure of the source are not as direct as they were before. The 
field interacts with itself due to the non-linearity of the theory, thus creating 
additional structure not attributable to the source. Examples of this are the 
critical phenomena in the scalar field collapse and the quasi-normal ringdown 
of perturbed black holes. Therefore, there is no direct relationship between the 
multipole moments of the source on the one hand and those of the gravitational 
radiation on the other. But irrespective of these issues, the null-datum in the 
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non-linear theory is ' 1'4 and it is the one that contains the radiative information 
of the gravitational field. The “news function” for full gravity has been defined 
by Bondi |^]) in a different way. It is related to the null-datum which is its 
M-derivative. In contrast to the null-datum which is a local quantity, the news 
is well defined only globally on > 7 + and therefore it is not as accessible as '^ 4 . 

6 Discussion and conlusion 

In all cases we have discussed, there is a null-datum which contains information 
about the structure of the field and its sources. These functions are those 
which can be registered on 7^+ by the infinite detectors. The extraction process 
consists of evaluating the news on 7'+ as a function of a Bondi parameter on all 
generators. Then one obtains the two polarizations of the radiation which arrives 
on We have described an idealised extraction process and we have argued 
that it is a natural idealisation of the real physical measurement. Furthermore, 
it is a procedure which is entirely intrinsic to a property which makes it 
very valuable within the framework of the conformal techniques in numerical 
relativity. 

Within these methods the extraction process becomes an unambiguous algo¬ 
rithm which allows the determination of the well defined news function up to the 
freedom of the choice of an origin of time for each detector and the indetermi¬ 
nacy due to the unknown boost between the source and the detector. The only 
approximation involved is the basic assumption that the arena be an asymp¬ 
totically flat space-time. In contrast to this, other extraction methods which 
are based on points in the finite region of space-time have to make additional 
assumptions in order to separate the far-held from the near-Held. 

We have shown in that one can accurately recover the news from T"'*’ so 
that the discussion of the present paper is not academic but a practical one. 
What remains to be tested numerically is the location of the generators of 7"’’’ 
on each time-slice and the determination of the Bondi parameters on them. This 
will be discussed in a later paper. 


References 

[1] H. Bondi, M. G. J. van der Burg, and A. W. K. Metzner, Gravita¬ 
tional waves in general relativity. VII Waves from axi-symmetric isolated 
systems, Proc. Roy. Soc. London A, 269 (1962), pp. 21-52. 

[2] J. Frauendiener, Numerical treatment of the hyperboloidal initial value 
problem for the vacuum Einstein equations. I. The conformal field equa¬ 
tions. to appear in Phys. Rev. D. 

[3] -, Numerical treatment of the hyperboloidal initial value problem for 

the vacuum Einstein equations. II. The evolution equations, to appear in 
Phys. Rev. D. 

[4] H. Friedrich, The asymptotic characteristic initial value problem for Ein¬ 
stein’s vacuum field equations as an initial value problem for a first-order 
quasilinear symmetric hyperbolic system, Proc. Roy. Soc. London A, 378 
(1981), pp. 401-421. 



REFERENCES 


13 


[5] S. Hawking and G. F. R. Ellis, The large scale structure of space-time, 
Cambridge University Press, Cambridge, 1973. 

[6] P. Hubner, Black hole spacetimes on grids with trivial boundaries, preprint 
AEI-062 jgr-qc/98040^ 

[7] J. D. Jackson, Classical Electrodynamics, John Wiley & Sons, New York, 
1975. 

[8] A. I. Janis and E. T. Newman, Structure of gravitational sources, J. 
Math. Phys., 6 (1965), pp. 902-914. 

[9] E. T. Newman and T. W. J. Unti, Behavior of asymptotically flat empty 
space, J. Math. Phys., 3 (1962), pp. 891-901. 

[10] R. Penrose, Zero rest-mass fields including gravitation: asymptotic be¬ 
haviour, Proc. Roy. Soc. London A, 284 (1965), pp. 159-203. 

[11] R. Penrose and W. Rindler, Spinors and Spacetime, vol. 2, Cambridge 
University Press, Cambridge, 1986. 

[12] A. SOMMERFELD, Vorlesungen iiber Theoretische Physik, Band III. Elek- 
trodynamik, Verlag Harri Deutsch, Thun, 4. ed., 1977. reprint. 





